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FRSFACE 


This  report  was  originally  submitted  as  a thesis  to  the  Faculty  of  the  Gradu- 
ate School  of  The  University  of  Texas  in  fulfillment  of  a requirement  for  the 
degree  of  Master  of  Science  in  Electrical  Engineering. 

Since  interest  in  the  cavity  resonator  was  occasioned  by  and  the  results  of 
the  work  are  pertinent  to  the  Office  of  Naval  Research  Contract  Nonr  375(01), 
it  was  felt  that  it  should  be  submitted  as  a technical  report  under  this  contract. 
Although  it  would  be  desirable  for  the  sake  of  completeness  to  extend  the 
calculations  to  larger  eccentricities,  the  main  interest  in  elliptical  cavities 
lies  in  the  small  eccentricity  region.  In  this  region  fall  the  eccentricities 
associated  with  unavoidable  deformations  of  cylindrical  cavities.  The  extreme 
teditxa  of  the  calcinations  dictates  the  use  of  approximations  in  the  calculations. 
The  particular  ones  developed  in  this  thesis  are  sufficiently  accurate  for  eccen- 
tricities less  than  0.4  and  are  subject  to  small  errors  for  eccentricities  between 
0.4  and  0.5. 


t 

i 

ABSTRACT 

Formulae  t .re  derived  for  the  quality  factor  ind  resonant  wave  length  of  an 
elliptical  res<  >nant  cavity  operating  in  the  TEnJ  mode.  Calculations  are  made 
and  curves  plotted  for  their  variation  with  change  in  eccentricity  for  values  of 
eccentricity  ljass  than  0.5.  The  necessary  integrations  are  numerical  using  sim- 
plifying asscjB?  tions. 

For  both  1 be  even  and  odd  inodes,  tne  quality  factor  increased  slightly  as 
the  eccentricity  was  increased  from  zero  to  a smaill  value.  A further  increase 
in  the  eccentricity  causes  the  quality  factor  for  t.ae  even  mode  to  decrease.  The 
range  of  eccentricity  (0  to  0.5)  used  was  not  sufficient  to  show  the  anticipated 
decrease  for .th'?  odd  mode.  The  eccentricity  range  considered  was  limited  by  the 
approximations  v-ad  in  the  method  of  evaluation,  l’he  approximations  were  con- 
sidered as  fully  Justified  for  eccentricities  less  than  0 .L  and  subject  to  seme 
. error  for  eccentricities  between  0.4  *nd  0,5. 
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CHAPTHt  I 


CnRODDCTION 

The  purpose  of  this  thesis  is  to  investigate  the 
effect  of  mall  amounts  of  elliptical  deformation  on  tho 
behavior  of  certain  characteristics  of  a resonant 
electromagnetic  circular  cylindrical  cavity.  Although  the 
circular  cylindrical  resonant  cavity  is  a special  case  of 
the  elliptic  cylindrical  resonant  cavity,  the  elliptical 
cavity  solution  cannot  be  expressed  in  terms  of  the 
cylindrical  functions.  The  Bessel  functions  used  for  the 
cylindrical  case  are  relatively  simple  and  numerous 
tabulations  are  available.  The  Hathlen  functions,  required 
for  the  elliptical  case  are,  however,  such  more  complex  and 
very  few  tabular  values  have  been  published.  The  resonant 
wavelength  and  quality  factor  in  the  elliptic  cylinder 
considered  as  a deformed  circular  cylinder  warrant 
Investigation  because  a physical  cylinder  nay  depart  from 
perfectly  circular  to  an  extent  determined  by  the  uanufacturinf 
tolerance;  external  forces  such  as  mounting  brackets  could 
also  cause  departure  from  circular.  It  seems  unlikely  that 
the  elliptical  cylinder  cavity  would  exhibit  such  decided 
superiority  over  the  olrcular  cavity  as  to  justify  the  considerable 
additional  manufacturing  difficulties  attendant  to  its  use. 

The  calculations  for  the  node  considered  here  show  no  advantages 
peculiar  to  tho  elliptical  modsa, 
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r.T». — iiw» 
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2 


s-  Calculations  ware  made  in  19U>  by  Ki Lnasr  and  Wiloon  ^ 

to  determine  the  variation  of  -wavelength  In  oertain  nodes 
with  the  elliptic ity  of  the  cylinder,  Rinser  and  Wilson  also 
derived  an  opression  for  the  quality  factor  far  ona  value 
of  eccentricity  far  the  TK^1  node;  this  thesis  will  coneidar 
the  TR111  node  with  several  values  of  eooentricity  for  both 
odd  and  even  excitation. 


1 J,  P.  linear  and  X»  0.  Wilson,  *Son  Basalts  on  Cylindrloal  . 
Resonators,"  Bell  Sfcrsten  Teehnical  Journal,  vol  26,  19U7,  p U2D 


CHAPTER  H 


(1) 


(2) 
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i 

ELLIPTICAL  COORDINATES  AMD  ELLIPTICAL  HAVE  PUMCTIOIB 

| 

(1)  The  Elliptical  Cylinder 

The  diaenal an s of  tha  cross  section  of  an  elliptical 

2 

cylinder  are  shewn  in  Figure  1*  The  quantities  2a  and  2b  are 
the  major  and  minor  axes  respectively)  the  fecal  distance  is  2q. 

The  perimeter  of  the  ellipse,  a,  will  be  kept  constant  when  the 
eccentricity  is  varied,  and  the  parameter  used  will  be  the  *averafa 
dimeter,*  D,  which  is  related  to  the  perimeter  by  the  formula i 

D s -£g^«g-  • JL. 

-TT*  ir- 

It  is  evident  that  for  tha  circular  ease,  D Is  the  diameter  of 
the  undlstorted  circle. 

The  eoaantricity,  <5  , is  defined  as  the  ratio  of  the 
aen&^fooal  distance,  q,  to  the  seoi-siajor  axis,  a.  The  eccentricity 
la  not  measurable  directly  and  there  are  two  other  directly  related 
quantities  which  are  often  used  Instead  of  the  eccentricity  as  a 

I 

measure  of  the  departure  frou  a circle.  One  quantity  is  the 
ell  Iptloity,  E,  ~efinedt 

E m between  major  and  minor  diameters 

major  diameter 

1 

a 

The  olliptioity  Is  related  to  the  eooentrioity  by  the  formula » 


These  and  all  other  symbolic  abbreviations  are  defined  on  page  . 11 

3 


7h®  other  quantity  •which  nay  be  used  to  express  the 
departure  from  circular  is  N,  defined  as  the  ratio  of  the 
minor  to  the  major  axis. 

N = b/e  *\il  - e2  , or,  after  an  elementary 
hyperbolic  trigonometric  identity  substitution, 

N a tanh  (are  seehQ) 

The  relations  (3)  and  (U)  are  plotted  in  Figure  3* 

Curves  for  the  variation  in  wavelength  and  in  the 
qualify  factor  are  plotted  against  the  eccentricity,  but  values 
of  either  S or  H can  be  found  by  using  Figure  3 In  conjunction 
with  the  curves  plotted  against  the  eccentricity. 

(2)  The  Elliptical  Coordinate  System 

The  elliptical  coordinate  system  is  ah  cam  in  Figure  2. 
The  orthogonal  coordinates  J and  q locate  a point  uniquely. 

The  elliptical  coordinates  arc  related  to  the  x-y  coordinates 
by  the  transformation  equations: 

x s q cosh  7 cos  ? 
y » q alnh  f sin  7 

The  equation  of  the  bound "-y  surface  of  the  ellipse  is: 
cosh  )T  0 c constant  • 1/  eccentricity 
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Tbs  eccentric ity , <?  , varies  between  sero  and  one  and  the  ellipse 
degenerates  as  6 — *1,  eoshF  — - if  f — *•  o to  a straight  lin* 
between  the  focii.  The  ellipse  degenerates  into  a circle  as  * — *0, 
eesh  7 «■  t — * «*•  The  elements  of  arc  length  in  elliptical 
cylinder  coordinates  are  da^  and  dSg  as  shown  in  Figure  3j  the 
s dlreotlon  is  into  the  paper  with  the  elaaent  of  arc  length  in 
that  direction  equal  to  one  da. 


(7) 


(8) 


If  the  definition  is  mades 

q^  c q (cosh2  T - cos2  7 
it  -will  follow  that 

da^  s q (cosh2?  - cos2?  )*  • q^  d 1 

ds2  a q (ooflftx 2JF  - cos2*;  )*  dy  d/j 

(3)  The  Pare  Equation  In  Elliptical  Coordinates 

The  tec  dimensional  rave  equation  In  elliptical 
coordinates  iat 

- k£  (cosh2?  - cos2 7 ) = 0 


This  equation  is  known  as  Uathien's  equation.  When  a product 
solution  is  assuasdf  and  the  Bernoulli  trial  method  of  separation 
is  followed,  the  equations  separate  to  give  two  ordinary  differential 


equations t 

2 

(ID)  (b  - q2  OOB2^  ) f ( 7 ) e 0 

2 

(n)  -*■  W - (b -k?  q2  ««h2I  ) f(f)  s 0 

®l  1 

where  b is  the  separation  constant*  The  solutions  of  the  first 
equation,  (10),  are  often  called  Mrthtsu  functions,  and  the 
velnt ions  of  the  second,  (11),  are  then  called  associated 
llathieu  functions.  Equation  (10)  transforms  to  (11)  under  the 
substitution  7 a ± i J , and  equation  (11)  transforms  to 
(10)  under  the  substitution  f c t i 7 , where  1 is  V-l. 

Solutions  exist  regardless  of  the  value  of  the  separation  constant, 
b,  but  the  solutions  are  periodic  only  for  certain  characteristic 
values  of  the  separation  constant • Some  authors  ^ consider  only 

the  periodic  solutions  of  Uathieu*s  equation  as  Uathieu  functions, 
but  more  recently,  Uathleu  functions  have  been  considered  as  all 
solutions  of  (9)  whether  or  not  the  conditions  for  periodicity 
are  satisfied*  ^ Uolstehlan  ^ has  an  extensive  discussion 
of  solutions  Where  no  restrictions  are  placed  on  the  separation 
constants.  In  the  calculations  which  are  required  in  this  work, 
only  solutions  which  are  periodic  in  7 will  satisfy  the  required 
boundary  conditions. 

3 

S*  T.  Whittaker  and  0.  N.  Watson,  A Course  of  Uodsrn  Analysis, 

Mtocdllan,  19U&,  p h0$ 

UTablss  Relating  to  Vathlsa  Functions,  The  Computation  Laboratory 
United  States  national  'duraau  of1  Standards,  1951 
5N.  W.  Uc Lachlan,  Theory  and  Application  of  Mathleu  Functions, 

Oxford,  19U7  * ~~  ' 

i 
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(it)  Matfaleu  Functions 

Msfchietx  functions  arise  in  stability  Irveoiigatlons  of 
various  mechanical  systerss,  tbs  theory  of  frequency  Emulation, 
loud-speaker  theory,  and  in  any  electromagnetic  or  vibration 
problem  •shich  must  be  stated  in  elliptical  coordinates* 

The  modern  theory  of  Matbieu  functions  i 3 credited  to 
Whittaker  and  much  of  the  aub sequent  theoretical  development  is 
credited  to  Luce,  Strutt,  and  Ma Lachlan.  A historically  complete 
list  of  226  references  is  given  in  !&Xaehlan«  ^ Because  both  the 
theory  and  numerical  ccnputaticns  of  Mathieu  functions  are  more 
difficult  than  those  fear  Bessel  and  Legendre  functions,  couplet* 
tables  chich  make  possible  the  actual  use  of  Mathieu  functions 
in  rsmsrUisJ  coaptation  have  lagged  far  behind  these  other  functions 
and  have  only  recently  become  available* 


The  elliptical  wave  guide  •was  first  investigated  by 

Q 

Cbu  in  1038  and  following  his  resrk,  the  first  numerical  tables 

o 

of  coefficients  far  Mathieu  functions  was  published*  Recently 
the  far  more  accurate  and  ccctensive  Tables  Relating  to  Math-leu 
Functions  have  been  published.  Notation  for  Mathieu  functions 


\ W,  MoLaehlsn,  Theory  gad  ApzdJ&atLcm  of  Msthiru  Functiona. 

^Ottferd,  l$h7  “ ' "** 

"SMSes  Rela^n^to Jfe|hi^^^^i^^K^C»o3yutatlQn  Laboratory, 

8L*  vbu*  "Electromagnetic  Waves  in  Elliptic  Metal  Pipe*,*  Journal 
of  Igpllad  Physios,  vel  9,  3$38»  p $83 
97*  A*  Stratton,  Ps  U*  Merse,  L*  J*  Shu,  and  H.  1*  Botnar, 
end  spheroidal  Ware  Functions,  Mar  lark,  Wlley| 
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Solutions  to  equation  (10)  can  bo  found  from  the  formula* 

^(0,00=7 ) : Id”  SSk, 

■where  SpQ(c,cos  9)  is  the  angular  hathiou  function,  p signifies 

either  • (even),  or  o (odd),  n is  the  order  of  the  function  (and 

for  the  node  considered  always  one),  c is  a parameter  which  is 

defined  later  ^ and  cos  9 is  the  argument  of  the  function* 

The  coefficients  for  the  right  hand  summation  are  found  in  Stratton 
12 

or  in  the  Tables.  The  cosine  functions  of  9 are  used  in  the 

suswatlon  for  the  even  functions  and  sine  functions  are  used  in  the 

sasasitan  for  the  odd  functions.  The  series  (12)  is  not  a Fourier 

series  because  the  coefficients  are  not  derived  from  the  Fourier 

13 

defining  integrals,  but  according  to  Stratton  they  apparently 
sst-ief*  the  conditions  of  convergence  necessary  far  term  by  term 
differentiation  or  integration. 


See  page  1? 

Xllcw#  cit.  p ?R  or  p 82.  Refer  to  Table  for  designation  used 
l^loc,  eifc.  tabulated  for  valuta  of  s.  Refer  to  Table  for  relation 
between  c and  s. 

^loc.  cit,  p 20 


table  is  reproduced  on  page  12.  With  the  exception  of  the  designation 
of  the  parameter  for  the  coefficients,  all  cf  the  notation  used 
here  agrees  with  the  notation  used  in  the  Tables.  Unfortunately 
it  was  impossible  to  use  the  Tables  for  numerical  calculations 
because  in  the  mode  chosen  for  investigation,  the  number  of 
values  of  the  parameter  in  the  desired  range  was  insufficient. 


- -----  :-r  .-xr — si . .ys  j 


, 1 


'i3»  enly  angular  functions  nhich  ■will  be  encountered 
in  ih©  rmxi*  will  bs  the  angular  function  ol  the  first 


kind  nith  solutions  of  the  form: 


S^(c,ooc^)  = ^ D*^  [(2*Ll)rl)  of  period  2 

c>£ 

oo^(c,oos  <7 ) s Eb^-,  ein  of  period  2 TT 

corresponding  radial  solutions  may  be  calculated  from 
a Joining  factor,  but  in  practice}  the  useful  expression  irtiieh 
converges  much  acre  rapidly  than  the  trigonometric  one  is  expressed 

✓ 

as  a sum  of  Bessel  functions: 

oo 

^ £ '-i>k  °*Li  ^ tM 1 > 

oO 

■fa^o.cosh  t ) ^ 'IV2  tanhl  X (o.CMh  f ) 


~7°\ 

order  m* 


is  the  Bessel  function  of  the  first  kind  and 


Eclating  v 

Statsiir&J? 


**3* 
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(13) 


(lfc> 


(1$) 

06) 


! 4*i 


The  only  angular  functions-  rbich  Trill  be  encountered 
in  the  T&iji  made  Trill  be  the  angular  faction  of  the  first 


kind  -sith  solutions  of  the  fora: 


S«Lj(c»oo8  rj ) - ^ ^®2kvl  S0S  of  period  2 TT 


So-,  (c,coa  *7 


> = 


sin  ifJikfrl)/?]  of  period  2 TT 


She  corr e sp ending  radial  solutions  may  be  calculated  from 
r Joining  factor*  but  in  practice*  the  useful  expression  Tdi.'-»h 
converges  much  more  rapidly  than  the  trigonometric  one  is  expressed 
as  a man  of  Bessel  functions : 


00 


J^CcjOosh  I ) 2 ^2  (-l)k  Ds^+1  ^ gbtO.  (cjCoah  I 


oO 


«fo^(e»©oah  l ) ” 4vi~  tankf  2_  {~l)"(2k)BoJk+1  ^ %i»L  (c*coaSl  £ ) 


Ik 


■where  is  the  Bessel  function  of  the  first  kind  and 
SX 

, lli 
order  an 


^Tabloa  Salat lag  to  Msthleu  Functions*  The  Computation  Laboratory* 
Oniiad  2t?ter  Bureau  of  ^andarcsV  ^9?1#  P sae 


f£S«LE  GT  WEATION  AND  CONVERSION  FACTORS 


ifctation 

Notation 

Used 

Used  in 

Here 

Tables  15 

2 

C 

s 

86^(0,0037) 

Ser(s,7) 

So?(c,cos7) 

So^Ca,^) 

K 

H 

K 

«fe«(c,eoohy) 

r 

Jer(s,?) 

Jo.  (cjQoshF) 

r 

Jo  (s,?) 
r 

be 

be 

r 

y 

* 

bo 

r 

bo 

r 

Rotation 
Used  Ixl  , 


3oi(c,cos^) 
r ' 

SOy(c,cos 7) 
* 


Jer(c,eosh|) 

J3^(c,coshf) 


Notation  Notation 

Used  in  17  Used  in 

liiLachlaii  ' Tanfp® 

l*q  Q2 

=er(f,q)/A  SeJ(c,eosy) 

8s_(7#q)/B  80^(0,0037) 

a£/a  I? 

4*  s 

Ce  (f,q)/Ag  (s)  Re  (c,coshf) 
r er  * 

Ser(?jq)/3gar(s)  P.or(c,coshf) 

a + 2q  b 

r r 

b <■  2q  b ' 

x*  r 


. loc,  eit*  p xxxviii  except  for  last  column 
U loe.  oit. 
loo^  oit* 

18  Cf  Co  Tang . ’’Propagation  of  ElectrcBaagnetic  Wares  in  Hollow  Metal 
Pipes  of  Elliptical  Cross-Section^"  19h9,  University  of  Texas  Thesis 


. -vt  ■'* 
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CHAPTER  HI 

a summary  of  orm  wcrk  os  elliptic  guides  ai©  cavities 

(1)  The  Results  of  Chu  and  T&sg 

The  problem  of  the  propagation  of  eleetrooiagnetic  waves 

in  hollow  pipes  of  elliptic  cross  section  has  been 

investigated  theoretically  by  Chu.  Chu  studied  the  six 

lowest  order  waves.  With  the  exception  ef  nodes  in  the  cylindrical 

nine  which  exhibit  circular  symmetry  (HI  , and  TE__),  when  the 

01  ©1 

cylinder  is  deformed  to  an  ellipse  both  even  and  odd  elliptical 
modes  are  generated  with  their  relative  magnitude  depending  on 
the  -polarisation  of  the  excitation*  Because  of  this  splitting, 
slight  deformation  of  the  cylindrical  guide  nay,,  unlike  deformation 
of  the  rectangular  guide,  lead  to  instability*  The  cavity  considered 
in  this  thesis  may  be  regarded  as  a very  short  wave  guide  shorted  at 
the  ends  so  that  the  generation  of  two  modes  does  not  lead  to 
instability,  but,  rather,  to  a broadening  of  the  frequency  response 
of  the  cavity  due  to  the  splitting* 

Ghu's  article  covered  the  theory  of  elliptical  wave  guides 
but  omitted  much  of  the  numerical  calculations  which  were  used  in 
obtaining  his  results*  These  numerical  calculations  were  reworked  in 
detail  by  Tang  in  1 9h9*  Both  Tang  and  Chu  obtained  curves  for  the 


I? 


L.  J.  Chu,  ^Electromagnetic  Waves  in  Hollow  Elliptic  Pipes  of 
Metal, " Journal  of  Ann  lied  Physics,  vol  9,  September,  1938 
20  c.  C.  Tar^'^Td^^tTO?ry  TSaSBTWteagnetic  Waves  in  Hollow  Metal 
Pipes  of  Elliptical  Cross-Section,"  19h9,  University  of  Texas  Thesis 
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A*'-  A" 


WBF&i 


variation  in  the  cutoff  wavelength  ms  a function  of  the 


eccentricity  for  a guide  of  constant  periphery.  The  curves 

1 * 1 1 1 M % 1 1 1 . 21 


21 

obtained  by  Chn  are  reproduced  by  Sarb&cher  and  Edson  and 
22 

by  Uoreno,  Both  Tang  and  Chu  obtained  curves  plotted  against 
eccentricity  far  roots  of  the  equations: 


Jo  (rei^)  a 0 required  fax  even  TU  modes 


Jo  (rol^)  s 0 required  far  odd  TU  modes 


& &Vn)  m 0 required  for 


even  TE  modes 


Jo  (ro^j[)  s 0 required  far  odd  TE  modes 


Unfortunately,  the  accuracy  required  in  the  numerical  work 


of  Tang  and  Chu  is  not  sufficient  fox  the  calculations  of  the 


resonant  wavelength  and  quality  factor  of  a cavity. 


(2)  Einaer  and  Wilson  Results  on  Cylindrical  Cavities 


Kinaer  and  Wilson  determined  the  root  values  of  the 


applicable  equations  (17)  or  (15)  correct  to  five  significant 


figures  for  nine  modes  in  the  elliptic  cylinder » the  even  TE  , 

Oln 


the  even  and  odd  TU^^,  the  even  and  odd  the  even  and  odd 

T*32n*  rran  ^oin*  even  T&,^n,  The  first  subscript 

Indicates  the  number  of  variations  in  the  angular  direction,  the 


ond  subscript  indicates  the  variations  in  the  radial  direction. 


Far  a resonant  cavity  the  third  sub  script  indicates  the  variations 


In  the  axial  direction,  but  this  does  not  affect  the  value  of  the  searas 


R.  X,  Serbaohsr  and  V,  A,  Edson,  turner  and  ultrahigh  Frequency 

SgiaMrirg,  Wiley,  19U3  

TTiCreno , Uicrcwgve  Transmission  Design  Data,  UcQraw-Hill,  19uo 


J,  F.  Eirier  and  I,  G,  Wilson,  “Some  Results  on  Cylindrical 
Resonators,"  Bell  System  Technical  Journal,  vol  26,  19U7,  p UlO 


I?  t 

: 


i 


«U.i  * 


\ 


r* 


Klnser  and  Wilson  determined  an  empirical  aquation  for 


the  ratio  of  tho  porimoter  to  the  cutoff  wavelength  for  three 


nodes  (tho  even  TEq^,  the  even  and  odd  Tit^)  as  a 


function  of  the  elliptlcity,  B* 


Xlnaer  and  Wiloon  derived  an  expression  for  Q for  one 


mode  and  one  value  of  eccentricity  (even  TEq-^  with  an  eccentricity 


of  0.U81U).  The  circular  synaetry  of  the  TE^n  makes  the  calculations 


necessary  to  obtain  the  quality  factor  simpler  than  those  necessary 


to  find  Q far  the  TByn  mode.  Their  article  does  not  give  any 


details  of  tho  methods  used  to  make  calculations,  but  the  sparse 


outline  of  method  of  calculations  indicates  that  the  procedure 


was  the  same  as  that  used  in  this  work,  with  the  exception  of  a 


different  formula  for  numerical  integration*  ^ Since  their 


only  result  for  Q is  one  point  on  a curve  for  a different  mode  than 


those  considered  in  this  thesis,  no  numerical  comparison  can  be 


made  of  results* 


See  juyje  36 
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OIAPTHl  IV 


(19) 


(20) 

(21) 


(22) 


DETBR1CDIAT ION  0?  THE  RESOHALT  WAVEL35GTH 


(1)  Field  equations  for  tho  elliptical  cylinder  cavity 


The  equations  for  the  components  In  an  elliptical  pipe 


25  26 

in  the  TE^  mode  aret 


H - 


E = 


H * = 


B 9p,( c,coa^)  ^(ejeoshf) 
0 

*> 


ik 


B = -B 


’l  kl 


l 3px(^)  ^ (f  ) ' k3*) 


H.  * 


*1 

UO  f*- 


* * -8  -^2j—  Sp'(7  ) ( t ) •i<wt  * k3,) 

"l11! 


27 


where  is  the  propagation  constant,  B is  a complex  amplitude 

constant  which  depends  on  the  relative  even  and  odd  nodes  excited,  w 

is  the  angular  frequency,  /*■  is  the  permeability  of  the  dielectric 

in  the  guide,  and  the  primes  denote  either  ■ & ■ or  ^ 

af  An 


The  boundary  conditions  require  that  E(  { ) s 0 where 


f0  is  the  boundary*  This  Ixplies  that 


i 


» • u ? 2 

Jp  (c,coahfo)  z 0 and  with  the  definition  rp^  = Q «►**»/*♦ 

ff  *l(^U)S  0 


^ TE  (transverse  eleotric)  is  often  written  as  H (since  only  non-scro 
oonpenent  in  the  s direction  is  H^) 

26  J*  A*  Stratton,  ELectrcBagnetlc  Theory,  UcQraw-IIill,  19U1#  p 375 

27  k-,,  the  propagation  constant,  Is"  of'ton  written^ 

28  In  all  the  following  work,  the  subscripts  and  primes  are  CKltted 
from  the  r designation  for  root  since  the  only  roots  which  occur 


will  be  either  r^^  or 


16 


17 


Combining  the  waves  that  travel  in  the  positive  s direction 
with  those  which  travel  in  the  negative  s direction  and  aaklng 
the  indicated  trigonometric  substitutions)  the  field  equations 
for  the  TK111  mode  in  the  resonant  elliptical  cylinder  are 
obtained.  The  tine  function  e is  suppressed  and  the  equations 
are* 

2 

Hz  s -B  3p1(c,cos»j)  J^(c,coahf)  sin(k^») 


H 


I * 


H ^ a 


E 


B!  = 


-i  B k-, 

-i  B k, 



% 

B k 
*1 

- B k 


1 

Sp1(c,cosy)  (c,coshf)  oos(k^z) 
Sp^CjCos^)  Jp-jO^cosh?)  cos(k^s) 
^(0,003  if)  Jjp|  (c,ooshfO  sin  (k^s) 
8J>1(c,cos7)  Jp^CjOoshf)  sin  (ky») 


E - 0 

s 


The  radial  and  angular  coefficients  are  tabulated  in  Stratton 

30 


29 


for  given  values  of  the  parameter,  c. 


C a 2TTq/A( 


^ Stratton,  Horse,  Chu,  and  Hutnar,  Elliptic  Cylinder  and  Spheroidal 

Wave  Functions.  Wiley,  191*1  - 

30  The  parameter  o should  not  be  confused  with  the  e often  used  to  designate 
the  speed  of  light. 
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The  perimeter  of  an  ellipse,  s,  is  related  to  the  aemi^focal 
distance,  q,  and  the  eccentricity,  (?,  by  the  formula: 


nr- 

Jt\ 


Q^COS^fJ  dq 


= Js-  e ((?)32 

e 

•where  E ( <2 ) is  the  complete  elliptic  integral  tabulated  in 
33 

Peirce  for  values  of  arc  sin 

It  irill  be  convenient  to  use  the  parameter  ^/s  used  by 
Tang  and  Chu,  and  c may  be  expressed  in  terns  of  that  parameter  by 
substituting  from  equation  (31)  into  (29). 


Vs 


(2)  Derivation  of  an  expression  for  resonant  wavelength 


The  parameter  used  to  express  the  shape  of  the  cavity  will  be 
R,  defined  as  L/d  where  L and  D are  the  length  and  "average  diameter" 
respectively  of  the  cavity.  The  following  definitions  are  made* 


k,  the  wave  number,  k ■ 27T/A 


s rS/q 
k2  = ^ + 
P = TT/2 


J This  s should  not  be  confused  with  the  s used  in  Tables  of  Ilathiau 
Functions  which  is  equal  to  c^,  nor  with  the  s used  in  ’Klnser  and 
Wilson  which  is  the  reciprocal  of  the  Aq/s  used  here. 

32  This  E ( $ ) should  not  be  confused  with  the  E used  by  Kinder  and 
Wilson  to  denote  ellipticity. 

33  3.  0«  Peirce,  A Short  Table  of  Integrals , Ginn,  1929,  p 121 
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It  follcsre  from  equations  (35)  and  (22)  that 
= (r/q)2  - ^2/<-c 

and  for  propagation,  k,  must  be  a pure  imaginary  and  '/*■ « 

2 

is  greater  than  the  quantity  (r/q)  , 

kj  : ik3  : i - (r/q)2 

- (r/q)2 

At  the  resonant  frequency,  ia  aero  so  that 

p 2 

(r/q)^  and,  after  solving  for 
Xe  - 2^Vr 
For  the  TE111  soda, 

k3  r iyt 

Hhen  this  value  far  is  set  equal  to  the  value  far  k^  from 
equation  (39),  and  the  equation  ia  solved  for  A,  the  resulting 
equation  1st 


'(1/2L)*  ♦ (r^tq)1 


When  the  value  free  equation  (Ul)  is  substituted  into  equation 
(k3)  and  both  sides  of  the  equation  are  divided  by  D to  make 


f 


the  resulting  expression  dimensionless,  the  equation  obtained 


(WO 


- ' '**•  - * ' - — a-  H 

» 
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is; 


The  values  of  the  parameter  X^/s  are  platted  as  a function 
of  eccentricity  in  Chu.  ^ It  has  already  been  noted  that  the 
values  of  \/s  obtained  by  Chu  are  not  sufficiently  accurate  to  be 
used  in  determining  values  for  the  quality  factor,  but  the  values 
for  the  resonant  wavelength  from  equation  (iii*)  can  be  evaluated 
very  easily,  without  the  use  of  Llathieu  function  tables,  if  Ghu's 
values  f>or  ^/s  are  used.  Figure  1*  is  a plot  of  equation  (UU ) 
using  Chu’s  values.  It  is  noted  that  these  resonant  wavelength  values 
are  not  sufficiently  accurate  and  do  not  enter  directly  into  the 

■v£  > 

values  for  the  quality  factor. 


If  a more  precise  determination  of  the  resonant  wavelength 
is  desired,  it  is  necessary  to  combine  equations  (31)  and  (1*1)  to 
obtain  the  relation  which  was  used  in  evaluating  the  quality  factor: 


(US) 


7 


(r)  ( v»>  » (V/2>  / E <e ) 


whan  the  eccentricity  is  equal  to  xero,  E (<2. ) is  equal  to 
TT/2  so  that  ^j/s  becomes  the  reciprocal  of  the  root.  ^ If  both 
sides  of  equation  (UU)  are  multiplied  by  D and  the  substitution  is 


??  loc  cit 

■*°  It  sh<*ild  be  observed  that  the  variation  in  resonant  wavelength 
dees  affect  Q since  it  appears  in  the  term  In  the  numerator  of 
equation  (82) 

37  Kinagr  refers  to  his  parameter  "s"  which  is  the  reciprocal  of  >»c/s 
as  the  "root  value  adjusted  to  the  eccentricity." 


i mi  i il» itai  mu  lirnMiiiii'illMiiiiiMWrtminni  M» 


c 
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(U6) 


Bade  for  zero  eccentricity  Jh*oo  equation  (1*.$)*  the  aquation 
obtained  1st 


2 

Ao1”  ' Vfr/fcS)*  * 0/L)5 


This  agrees  with  the  fcsnsula  far  wavelength  foe*  the  circular 

38 

resonant  cavity  given  by  liontgoiaery. 
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C«  Q*  Uontgcraery,  ed«.  Technique  of  Mlnrowaww  Meagre 
Radiation  Laboratory  ferSsI,'  IctrS^SSli,  i^'f % 2f7 
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chapter  v 


IBPmilNATIOlI  of  the  quality  factor,  q 


(l)  Derivation  of  an  aanregglan  for  the  anal  Ity  factor.  ^ 


The  quality  factor,  Q,  of  a resonant  cavity  is  ordinarily 

defined  as  the  ratio  of  the  product  of  the  angular  frequency  and 

39 

the  energy  stored  to  the  average  power  loss,  The  total  energy 
in  the  electric  and  magnetic  fields  remains  constant  (neglecting 
losses)  and  the  maximum  electric  energy  is  equal  to  the  maximum 
magnetic  energy.  It  is  sufficient  to  consider  either  the  electric 
or  magnetic  stored  energy.  Since  one  component  of  the  electric 
field  is  absent  in  the  TE111  i node,  it  is  moire  convenient  to 
consider  the  stored  electric  energy.  This  can  be  found  by 
integrating  the  E j and  E ^ components  of  energy  over  this  volume. 


Q = 


The  stored  electric  energy. 


( f\z if  /£= 

/«J*0  Jfm 0 J sao 


W !tj  ds^ds^ds^ 


The  differential  element  in  the  z direction,  ds^,,  is  equal  to 
da,  and  the  only  function  cf  z nrhich  appears  in  the  expressions  for 
E ^ end  2^  ^0  is  sin  ('k^s)  or  sin  (tra/L)  so  the  integral 


^ 3,  B*mo  and.  J,  R.  Fhinaery*  Fields  and  Waves  in  Modern  Radio, 
Filey,  19iuj,  p 378 
^ Equations  (25)  and  (26) 
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21* 


with  respect  to  b is* 


n.o> 

VH/  / 


(ire/L)  ds  a 1/2 


This  reduces  the  expression  for  stored  electric  energy  to  the 
double  integral: 


(50) 


I v**  /K>  _L* 

Jfd)  J$aO  u 


Substituting  f?oss  equations  (8)  and  equations  (26)  and  (2?)* 


(5D 


Uc  3 


Clv°  ([spi(  7 )]  2 &( E ^ 


2te 


7«2/rrfBr  LeB^kVr  . 


da. 


/T*6"  f 

rL  A 


ysO  -/^aO 


U€ 


The  radial  and  angular  functions  are  entire  functions  of 
^ and  y respectively,  anH  <*ach  double  integral  nay  be  written 
as  the  product  of  two  single  integrals.  It  is  noted.  fro» 
equation  (8)  that  do,  is  equal  to  q1  &?  and  is  equal 
to  d f . tfhen  these  substitutions  are  made  and  the  terras 
indspendest  of  the  variables  of  integration  are  reasored,  the 


« 


pjWMMM  ^fcaWWBWWIfiNiiW^  1 * 
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expression  far  the  stored  energy  can  be  -Brittens 


x-  B^k2  n 32 


R«aa.i 


*LT^r‘t 


Ths  fcllvsing  abbreviations  for  integrals  Trill  be  adept  edt 


/7«2TT  « 

^ * L*  C*l(  7 5 2 d? 

1*1*0 


^ = /IT 


“>  = dI 


nip  = rr*^[^  (?  )js  dT 

•* 


These  integrals  cannot  be  evaluated  analytically  and  the  evaluations 


aeurfc  be  made  by  a cccfeination  of  integration  of  series  and  numerical 


methods*  The  actual  evaluation  of  the  integrals  Is  discussed  in 


Chapter  VI* 


The  expression  for  the  stored  electric  energy  is  written 


using  these  abbreviations  ass 


L B^k2 


[ipUp  ipnjp] 


■ -4-W 


•■  i ' — L‘  i 


I ME 


- 


(58) 


159) 


(60) 

(61) 


(62) 


The  power  loss  in  the  resonant  cylinder  is  doe  to  the 
copper  losses  from  the  currents  flowing  in  the  side  -sell  and  in 
the  end  plates.  The  average  power  loss  is* 


/. 


'surface 


V 


•where  R is  the  surface  resistivity  as  defined  in  Rama  and 

s 

Whinnery.  In  the  side  wall  the  I component  of  current  is 
sere,  since,  by  equation  (2k), 

Hr(S0)  = 0 

The  other  two  components  of  current  present  in  the  side  wall  are 


evaluated  with  ¥ « and  when  squared  are* 


U2 


Hg~  - B2  ^{^(7  j]2  2 sin2(k3«) 

v = [^i<7  >j2  r*jM2  °«2<k3*) 


These  values  must  be  integrated  over  the  side  wall  which  requires 
integration  from  a * 0 to  8 » L and  from  <j7»  0 to  *jm  2^i 

Tj 

The  value  of  the  integral  in  equation  (U9)  is  the  same  for  an 
integrand  of  either  the  sine  or  cosine  function  squared  so  that 
the  integration  over  the  s range  of  either  of  equations  (60)  or  (61) 
yields  1/2  * 

The  integral  of  Hg  aver  the  side  wall  may  now  be  written* 


]Wi 


B2 


s 

side  wall 


f2tr 

Jo  [*a <i>J' 


lta* 


?y***"—  « mmmmm  ■ —— *m  m 

loc*  cit»  p 209 

ll2  from  equations  (23)  and  (2£) 


I 


nSStJ*  « »<W  . 


The  integrand  is  multiplied  and  divided  by  q.,,  and  the  substitutions 
from  equations  (7)  and  (8)  evaluated  at  F*  f Q 

x q(coah2?0  - eos2^)^ 
d 7 S ds2  / qx 
are  made  so  that  (63)  becomes  t 


sldssrsil  2 

Factoring  out  q cosh|Q  from  the  integrand  and  substituting  Q s l/ooah?o 
in  the  integrand,  equation  (65)  becomes 

sidewall  2 70  ^ L 1 /J 

| a 

A further  abbreviation  is  made* 

ir^  s J (i  -SSeoe2^)^  [35^(7 )]2  d^ 


and  then  the  integral  of  H over  the  side  trail  is  written* 

% 

J I iy  2 = i B2  [^iOfojl2  1 9 ““*>?<>  ™P 

sidewall 

The  integral  of  the  square  of  the  7 consonant  of  H over  the 
side  wall  is  written* 


sidewall  2 -/0  ^ 1 q^ 


Admittedly  there  is  a plethora  of  integral  abbreviation  symbols. 
Hcwsver,  the  final  expression  for  Q would  require  three  pages 
If  written  without  these  abbreviations,  so  they  must  be  accepted 
as  a necessity  rather  than  a confusing  convenience* 


9«M&&5*e  Hs&txae 


- - 
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The  substitutions  from  equations  (63)  and  (6U)  are  made  to  gives 


sidewall 


which  becomes 


0 q(cosh  | - cos  n) 


/lH,l 

** t i 


sidewall 


B"  ^ fe  Cf0)J‘ 

2q  coshr 


-jf  (i  a[^i^7  3 2 d7 


when  l/q  cosh  f .La  factored  out  of  the  integrand  and  the 
o 

substitution  $ ~ l/cosh  'tQ  is  nade  in  the  integrand. 

The  further  integral  abbreviation  is  mads* 

IVS>  ajo  (1  -SW^)  [Sp-J^)] 

Then  the  integral  of  Hrj  can  be  written! 

f 1h,|2  = f-lfeiJyJ.2 1 j4 

3iderall  2q  cosh  YQ 

Substituting  from  equations  (68)  and  (73)  into  the  expression 
for  average  powe"  loss,  equation  (£8),  the  expression  for  average 
power  loss  in  the  side  wall  is  obtained* 

P1OT  » R«  t^A?)2  [*4 5 CS3hfo  =* 

This  may  be  put  in  a somewhat  more  convenient  forms 


It  is  seen  from  equation  (23)  that  the  s' consonant  of 
H is  zero  at  the  end  walls  where  z « 0 or  8 »7T«  \7hen 
equations  (21;)  and  (25)  are  evaluated  at  either  of  the  end  walls 
the  square  of  the  current  in  one  wall  is? 


lK,\2  = ~ 2-  [sp^T)]2  [^(nj2 


When  the  substitutions  ds^  r dl  and  ds^  = q^  dy  are 
made,  the  integral  equations  become? 


■ 0V7>J2  L*i(r)i2d7dr 

The  integral  abbreviations  stated  in  equations  (53)*  (5U)* 

(55)*  and  (56)  are  used.  The  fact  that  there  are  two  end  walls 
provides  a two  which  cancels  the  factor  of  one-half  in  equation  (58) 
so  that  the  total  power  loss  in  the  end  walls  becomes? 

PL  - B2  k?  R (Jpn£  + 3jpI3p) 

ew  3 s 


. ■■ 


iwHMWd 


liquations  (£7)#  (7U)»  and  (79)  provide  the  information 


required  to  substitute  in  the  equation: 


‘,jJ  U 


Average  Power  Loss 


The  angular  frequency,*-^  , o an  be  expressed  in  terms  of  k by  using 


the  relation: 


2 IT  _ k 
“ vi/Ts 


from  equation  (33). 


This  gives  the  formula  far  Qi 


^ i P (3j?n£  + ^nj>) 


Rs  _L 


It  ulll  be  convenient  to  considar  Q <5 /a  instead  of  Q*  The 


relation  used  -will  be* 


**  ^ = — A,  then  Q <5/>>  a A/tT 


Also  the  substitutions  are  made  from  equations  (36)#  0*2)#  and  (3U)t 


p txTT/2 


1/2  s pA3 

q/e  s rAx 


This  trill  make  diroot  comparison  -with  published  curves  far  Q 
for  the  circular  case  possible* 

Hamo  and  Whinnory#  Fields  and  Waves  in  Modem  Radio#  Wiley,  19l*h,  p 211 


. .'5-'a4fs3a*t #&■«««*-  ‘ **■*&■* ' 


,iK-  .*%»'  ."WSW*#1-!  *H***Miil!X 


(82) 


1 

2 tt 


n-» 


Q 


k3  (ipni  + 3^np) 


— J 2[r2UIp  + ( if)17? 

2?  i_  \ j_ ' 


p- — — - 

+ — (Spnp  + ipn?) 


Another  integral  abbreviat ion.  is  madei 


(83) 


7p  s 3j>I3p  + 3pl3p 


(8U) 


From  equations  (33)  and  (UU) * 


k a 2-TT/x  a ~ 


(p2*2  + 


£3] 


* 


and  from  the  equations  on  the  preceding  page  and  the  relation  R a D/L 


(85) 


i -- 


8 vh.3 


and*  finally,  using  equations  (3U)  and  (33) 


(86) 


a r e/q  ■ Ip*  E(  e )/s  a 


D V3 


Whan  the  substitutions  from  equations  (83),  (3U),  (85)#  and  (86)  are 
made  in  equation  (82),  and  the  numerator  and  denominator  are 
simplified,  the  equation  for  Q used  for  making  calculations 
is  obtained. 
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CHAPTER  VI 

CALCUUTiCfN  OF  THE  QUALITY  FACTOR,  Q 
(1)  Detgnalnation  of  the  roots  tad  gagnfc&SflS  *hlch  fniinw  dlrftntlv. 

The  first  step  in  the  calculation  of  the  quality  factor,  Q, 
ie  the  determination  of  the  roots  of  the  boundary  condition 
equation  (22)  far  the  values  of  c which  ere  used*  It  develops 
that  far  snail  values  of  eccentricity  (between  aero  and  0.5)  that 
the  value  of  the  parameter  c varies  between  *ero  and  one  in  this 
node*  The  new  and  more  accurate  Tables  Relating  to  llathieu  Functions 
cannot  be  used  because  the  only  values  of  coefficients  falling  in 
the  desired  range  ere  for  values  of  c of  xero,  0*707,  and  1*0.^ 

The  Elliptic  and  Spheroidal  Ware  Functions  provide  the 
coefficients  far  values  of  o at  intervals  of  0.2  accurate  to  five 
significant  figures. 

ho 

Combi' s Eleven  Place  Tables  of  Bessel  Functions 
were  used  to  evaluate  the  Bessel  functions.  The  coefficients  for 
the  llathieu  functions  are  sero  beyond  far  the  range  of  c 
considered  in  both  of  the  series  (15)  and  (16)  so  that  Bessel 
functions  of  the  first,  third,  and  fifth  order  were  the  only 
ones  required.  It  is  desirable  to  determine  the  root  values  to 
five  significant  figures,  but  the  arguments  far  Bessel  functions 


h 6 


50 


tar 


lco.  oit.  pages  58  and  155 


{*7  values  of  s of  aero,  0.5,  and  1.  Refer  to  table  on  page  12 
h 0 loo.  cU. 


k?  E.  Combi.  Eleven  and  Fifteen-Place  Tables  of  Bessel  Functions,  Dover,  191*8 
50  The  same  interval  is  used  in  Jahnke  aixf~5a5e,~Tables  of  functions, 

Dover,  191*5 
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are  given  in  Csuabi  to  increments  of  0*01  in  the  neighborhood 
of  the  roots#  To  determine  the  roots,  auxilary  tables  were 
made  using  Linear  interpolation  for  Bessel  functions  of  the 
first,  third,  and  fifth  orders  in  increments  of  0*001  from  1*80 
to  2*10.  This  gave  values  of  the  needed  Bes3el  functions  for 
three  hundred  rallies  of  the  argument  in  the  desired  range* 


Use  roots  of  the  radial  functions  are  given  to  an  accuracy 
of  0.01  in  Tang.  For  the  even  mode,  using  the  roots  of  Tung 
as  a starting  place,  values  of  the  radial  function, 

oO 

-Je  (c  cosh?)  a \^/2  2 (~l)k  Dn^i  if  2k+l  Cosh1£) 
x k«0  'J 

were  plotted  against  c cosh  X and  a maximum  was  found  graphically. 

This  maximum  occurs  at  the  root,  r,  (which  is  equal  to  c cosh  ?Q). 

When  the  root  is  found  for  a given  value  of  o,  then  cosh  Crt  and 

oar  be  determined,  and  the  eccentricity  is  the  reciprocal  of  oosh  ? o- 

For  the  odd  mode  the  roots  must  be  determined  saaewhat 
differently  since  the  term  tanh  ? appears  in  the  formula  for  the 
radial  function, 

Jo^Os  cosh?)  = fry!’  tanh^  ( ~1 )k ( 2k ) (c  CO^i?) 


This  function  was  plotted  against  values  of  • and  the  maxinaan  determined# 


The  V?FA  Tables  of  Circular  and 


Sines  and  Cosines  was 


£7  loc.  cit*  page  !i3 
52  Tables  of  Circular  and  Hyj 
AfcnSstratibn,  -lew  YorkT 


>lic  Sines  and  Cosines,  Work  Projects 


3£ 


V9»d  to  evaluate  Jf.  Stops  of  the  argument  are  0.0001  so 
that  no  interpolation  was  required*  When  Tn  is  determined, 
the  eocentricity  and  the  root,  r,  can  bo  determined  by  a process 
which  is  the  reverse  of  that  used  for  the  even  mode* 

Kinasr  and  Wilson  ^ give  root  values  far  the  oven  node 
to  five  significant  figures.  They  make  no  conmant  on  the 
probable  error  in  these  figures;  this  would  lead  to  the  assunptioD 
that  a more  accurate  Interpolation  formula  vas  used  than  a linear 
one.  The  even  roots  obtained  by  this  author  by  the  method  outlined 
above  agree  with  the  Klnser  and  Wilson  values  within  one  ten-thousandth. 
■Actual  calculations  were  made  using  tho  Klnaer  and  Wilson  values. 

Klnaar  and  Wilson  did  not  detormina  values  for  the  odd  mode  of 
the  ESod*«  According  to  Scarborough  ^ these  root  values 

obtained  using  linear  interpolation  should  be  considered  as  subject 
to  a possible  error  of  two  ten-thousandths*  It  will  be  soon  later  ^ 
that  the  error  present  in  the  root  value  dominates  all  the  other 
errors  present  in  the  values  far  Q* 

The  values  of  the  ooqplete  elliptic  integral  from  equation 
(13)  are  tabulated  In  Peirce  ^ with  arc  sin  C?  as  tho  argument, 
linear  interpolation  was  used  to  find  the  value  of  E (0  with  a 
possible  error  of  0*0000£* 

gT’  ' ' 1 

\r  leo*  oit.  page  U28 

Soarbevugijt  Httwrloal  Analysis,  LSeOraw-Kill 

I'  Chapter  VII 
5®  loc.  cit.  page  121 
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I text  the  parameter  Xg/a  vu  determined  using  the  formula 
derived  from  equation  (u5) 


r E e ' 


All  of  the  roots  and  the  quantities  which  are  found  directly 
from  the  roots  are  tabulated  at  the  end  of  this  chapter. 


(2)  Evaluation  of  the  integrals 

(2TT  2 

1)  Evaluation  of  Ip  * I [3p(  rj  >J  d7 

-/n 


From  equation  (13): 


Se1(c,cos7)  I De^.^  cos[(2k*l)'|j 


This  aoqpreeslon  can  be  squared  and  integrated  term  by  term,  Because 
of  the  orthogonality  of  the  trigonometric  functions,  cross  products 
will  give  i iro  over  the  full  range  and. 


cos^  ®c  dx 


-L 


2 TT”  2 
sin 


dr  sir 


9 


so  that  the  value  of  the  integral  will  be  the  sum  of  the  squares 
of  the  coefficients.  That  is 

X.  ; TT  [(D^)2  - (D.1)2  - (D.|)2  - • * '] 

$8 

The  De  are  tabulated  in  Stratton,  horse,  Chu,  and  Hutner. 
r 


^ B.  ft.  P«lre«.  A 3hort  Table  at  fornula  U89 

5*  loo.  clt.  page~78 
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It  ahodd  bo  noted  that  far  calculations  In  which  it 

is  possible  to  use  the  Tables  Relating  to  Uathlau  Functions^ 

this  value  can  be  found  directly  from  tho  tables  tabulated  as 

N r tt/a2. 

r 

The  Integral  far  the  odd  function  differs  only  in  the 
trigonometric  function  involved  and  the  numerical  value  of 
the  coefficients*  The  Dor  are  tabulated  in  Stratton  ^ on 
page  82* 

(91)  Zo  Z ^ [(Dol)2  ♦ (Do3)2  ♦ (Do$)2 


These  integrals  are  also  evaluated  directly  in  the  Tables  Relating 

| p 

to  Uathieu  Functions  for  values  of  c as  N_  = Tl/B  . 

— - ■■  ”■  r 

t /2  ^ | 2 

ii)  Evaluation  of  3p  * J £sp^(  7 )J  d 7 

The  expression  in  equation  (13)  can  be  differontiated  term  by 
term  to  gives 

oo 

(92)  3e^(c  cos  7)  3 X.  -Oe^^  ^ 2k*1 ) sinf(2I*0.)7^ 

for  the  even  node,  and  for  the  odd  modes 

o o 

(93)  (o  0007)  ss  (2kel)  ooe[(2ksJ.)7j 

The  De^  in  equation  (92)  and  the  DOy  in  equation  (93)  are  the 
same  coefficients  as  thoee  in  equations  (90)  and  (91)  respectively* 
The  values  of  the  integrals  are  therefore! 


PT  

loo*  oit*  as  ?2*  See  table  on  page  12 
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(9ii) 

(95) 


(15) 


(96) 


I « = TT^Do^)2  + (3De^)2  4-  (5Db^)2  + * * *] 
ii  = TT[(Do[)2  f (3Do^)2  ♦ (5Do^)2  * * * •] 

ill)  Evaluation  of  lip  = j^°  (c  coshfj2  d f 

Proa  aquation  (15) 

oO 

Je,  (o  ccshf)  r \fiy2  2L.  (-l)k  De^^  $ (c  cosh?) 


For  a given  a,  the  interval  to  the  corresponding  jT  was 
divided  into  2h  parts,  and  values  for  Je^(c  cosh?)  were  calculated 
far  each  of  these  points*  Each  of  the  values  was  squared  and 
these  values  were  integrated  numerically  using  Sinpaon’s  one-third 


rulej 


61 


E.7EGRAL  = 


4 U(yx  4 7j  + •••*♦  j2y  ■*2(^2  4 7h  4*“43r22^  4 72lt 


where  h is  the  difference  botvreen  successive  abscissas  (that  is 

?/2li  in  this  particular  case)  and  the  y's  are  the  ordinates*  The 
o 

combination  of  the  error  in  the  coefficients  and  the  Bessel  functions 
was  calculated  to  be  0*00005*  The  inherent  error  in  3 inpeon*  a form  la 
far  integration  is  given  asi 


where  f^T  is  the  fourth  derivative  of  the  function  and  ray  be 


sr 


L.  Scarborough,  "freer ical  Analysis,  page  176 


approximated  by  the  fourth  differences.  The  fourth  differences 
are  of  the  order  cf  C,OCOC2  so  that  the  inherent  error  is  clearly 
negligible  with  the  coefficients  used, 

Kinzer  and  Wilson^  used  Weddle’s  formula  for  integration. 

The  inherent  error  is  considerably  less  than  the  error  with  Simpson's 
formula  when  the  integrated  function  has  sixth  derivatives  that 
are  a great  deal  less  than  the  fourth  differences*  Since  the  value 
of  the  coefficients  are  accurate  to  only  five  significant  figures, 
the  additional  difficulty  in  using  Weddle's  rule  seems  unjustified. 

The  estimated  error  in  these  integrals  is  twice  the  sum  of 
the  product  of  the  root  and  the  error  in  determination  of  the 
functions,  and  idle  product  of  the  average  of  the  function  value  and 
the  error  in  finding  the  root.  This  error  was  calculated  as  0*0009, 

The  integral  lie  is  calculated  in  exactly  the  same  manner  as 
lie,  using  equation  (16) 

•Jo.,  (c  cosh  f)  - VT/?  tanh?  (~l)k  (2k)lk>^^Q  2k+l'C  ccsh^' 

ksO  V 

s'C  I 2 

iv)  Valuation  of  Up  s j ° jjJp^  ( If  )J  df 

The  expressions  (IS)  and  (16)  are  differentiated  and 
evaluated  at  each  of  the  2?  points  which  were  found  in  iii).  The 
functions  are  squared  and  the  same  integral  formula  is  used  to  perform 
the  integration*  The  error  calculated  was  0,000^  which  is  larger  than 
the  calculated  error  for  the  integrals  lip  because  of  the  error  in 


finding  the  deriratiTB. 

y)  Evaluation  of  Hip  a f (1  - e2coa2fjft  [9p(  j]2  d^ 

Jo 

The  tern  (1  - e^ooe2  ^).  may  be  expanded  by  a binomial 
expansion  to  giro 

£l  - (1/2)  S2cos2»J  - (3/8)^008^  - • • 

All  but  the  first  two  terns  of  the  expansion  are  neglected  and 
substituting  from  equation  (13) 

Ills  s (1  - ^eoe27  ) [Se(7)J2d7 

If  the  trigonometric  substitution 
cos27  s £(1  + cos  2?) 

Is  made,  then 

me  = f2ir/l-i2(l^oos27)](_3e(7)]2d7 

•'0 

o /2fT  2 ,2^  2 

in*  = (i  - $€*)  (s«( *7  d<J  “ ^ ^ 3 008  27  d7 


The  integral  expressed  by  squat.  ?.jn  (100)  is  evaluated  In 


IJcLachlan. 
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oo 

Hia  . (1  - i^)TT  - f 2 * ZtD^XD^)] 
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loc.  cit.  page  79 


1 


4 f 


,>  I 


(102) 


(103) 


OoL) 


(106) 


A relation  for  the  odd  node  is  obtained  in  the  sane  manner 
as  fair  the  even  model 


nio  = (l  - ie2)  j2  [ so(^)]2d7  “$*f0  l?°(  9 )J  2 cos  27 

The  integral  (102)  is  also  evaluated  by  He  Lachlan 

» U - i*2)"'  - je2[-«D^)2  * Z (o^Xn^) 
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oO 

z 

lao 


/2Tt  Jp,  2 

vi)  Evaluation  of  IVp  r J (l  -.02cos^)  * [^>(7  )]  d7 

As  in  ▼),  a binomial  expansion  is  used  -which  in  this  case 
gives 1 

1 f |^ooe2<7  + ooa2q  - • • • 

Ifeglecting  all  but  the  first  two  terns  of  the  binomial  expansion. 


(105)  IVe  i j (1  f ^C2)  [se(  *7)J2  dy  4^  £e2cos  27  £e(  7 )] 2 d? 


which  is  evaluated  by  Uc Lachlan  as 

IVe  s (1  f ie?)TT*  TTie2  j^(Do^)2  5^  De^1D82k+3(2k41)(2kf3) 


Similarly-  for  the  odd  oaae. 


OO 


(307)  I7o  a (!4i«2)lT  + VK[*(&>i)2  + Z DoLflD°Lf3( 21041  )(2k^3) 

*•  kaO 


ST 


2cc*  oit.  page  79 


When  the  eccentricity  of  the  ellipse  is  aero,  the  constant 
c is  also  zero.  All  of  the  coefficients  are  zero  except  the  first 
one  so  that  the  angular  functions  reduce  to  trigonometric  functions 
and  the  radial  functions  reduce  to  first  order  Bessel  functions* 
Reference  to  equations  (90)#  (91)#  (9U)#  and  (95)  shears  that  the 
integrals  Ie,  Io,  Ie,  and  Io  all  will  have  the  value  ^ at  More 
eccentricity*  Further,  consideration  of  equations  (99)#  (102),  (105) 
and  (107)  shoes  that  Hie,  IIIo,  IVe,  and  IVo  have  the  value  at 
zero  eccentricity* 

Substituting  from  equation  (83) 

Vp  m TT(Up  + Up) 

and  (Up  + Up)  is#  in  the  degenerate  case: 

+ [}ji<f)32  af 

Equation  (108)  is  a fecial  farm  of  a Lasmel  integral  and  may 
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be  evaluated  directly  to  give : 

Vp  . 4-  r [^(r)}2  [l-  (Vr)2] 

The  ralue s for  tho  integrals  are  substituted  into  equation  (87) 
and  the  substitution  for  the  circular  case  from  eepatien  (U5) 


1 1 1 


TABLE  I 


Tbe  Root  Value 3 and  The  Associated  Quantities  for  the  Even  Hode 


<2 

>^8 
Hie 
I Ve 

Ie 

li 

He 

ui 

C*(r)]2 

Q 
R2 
R3 


0 

1*6102 

0.0 

0.  $13121; 

3.11059 

3.11059 

3.11059 

3.11059 

0.44689 

0.191 ilO 

0.53102 

0.20U56 

0.214702 

9.61078 


0.2 

1.8106 

0.10860 

0.530Ua9U 

3.12758 

3.14621 

3.14976 

3.14979 

0.45821 

0.17233 

0.53331 

0.205813 

0.216978 

9.61;996 


0.4 

1.8100 

0.21704 

0.549164 

3.086U; 

3.15931; 

3.17320 

3.17383 

0.48853 

0.13689 

0.53819 

0.206381; 

0.221071 

9.616895 


0.6 

1.8452 

0.32516 

0.557018 

3.01715 

3.17923 

3.  21401; 

3.21737 

0.51456 

0.10291 

0.54616 

0.207941 

0.236928 

9.55305 


0.8 

1.8484 

0.43280 

0.568021 

2.92922 

3.20025 

3.27436 

3.28523 

0.54168 

0.06451 

0.55780 

0.209873 

0.254519 

9.48998 


1.0 

1.8527 

0.53975 

0.585972 

2.74550 

3.23751 

3.35585 

3.38372 

0.47966 

0.03499 

0.56722 

0.192042 

0.291053 

9.06385 


TABLE  II 


V 


/ 


The  Root  Valnee  and  The  Associated  Quantities  for  the  TEm  lied* 


0.1* 

L.9082 

0*2113 

0.53006 

3.121*82 

3.19U03 

3.23791 

3.23851 

0.39611* 

0.23591* 

0.52351 


<5 

v» 

IXEo 

t 

IVo 

lo 

Io 

no 

no 

[Jo(rfl  2 

ooef  Q 
coef  R2 
ooef  R^ 


0 

1.81*12 

0.0 

0.5U312 
3.1 1*159 
3.1 1*159 
3.11*159 
3.11*159 
0.1*1*689 
0.191*1*3 
0.53102 


0.2 

1.8751 

0.1063 

0.531*77 

3.13723 

3.151*85 

3.16550 

3.16552 

0.1*3561 

0.21650 

0.5287U 


0.208682 

0.291*606 

0.1*81811 


0.8 

2.0075 

0.1*000 

0.52001 

3.0095b 

3.32977 

3.51*321* 

3.551*11* 

0.31817 

0.32391 

0.1*9973 

0.231*021 

0.1781*33 

0.510163 


1.0 

2.0751 

0.1*809 

0.51355 

3.07606 

3.1*1290 

3.79099 

3.81876 

0.27810 

0.36353 

0.1*8019 

0.253380 

0.167672 

0.532039 


0.201*56  0.203011 
0.211*702  0.201811* 
9.61*278  0j*30000 


0.6 

1.9560 

0.3061* 

0.52399 

3.10826 

3.25186 

3.36122 

3.361*55 

0061*58 

0.29106 

0.51379 

0.218826 

0.185251 

0.50107U 


h$ 


L.2 

2.1651 

0.5531 

0.1*91*97 

3.07210 

3.1*9823 

lu  11913 

U.17998 

0.21*357 

03*1521* 

0.1*5900 

0.281*1*92 

0.153956 

0.531*851 


3g§gjlig£giiiii 

^a»E»B!ggiBiaa 


k ■»»•!<• 

■•••>  •••••• 

ItflkyMlMItf  > • * •»  . t f t-  7» t « ••••-  . ♦ ' ■ 

k*  •••  '%  i — xx 


pOlII 

t»n«i»n«  I 

annua:  I 
nh:::n: 


i®i$SSi3S!Hi^d 

sarasgisraiBB 


•«*B5E  u£»3Bt  trIH  •«•*  n*j.  aH\  «:n  ton 

aiy  can  nrr  r?p  icc  ujri 


TABUS  HI 

Table  of  Values  of  the  Quality  Factor,  Q <£/*,  for  the  IR^  Modes 


Bren  Mode 


0.0 

0.2 

0.1* 

0.6 

0.3 

1.0 

o.o 

0.20U56 

0.20581 

0.20638 

0.20791* 

0.20987 

0.19201* 

o.5 

0.23529 

0.23670 

0.23752 

0.23951* 

0.21*210 

0.22322 

1.0 

0.26968 

0.2707!* 

0.27101 

0.27191 

0.25509 

0.25509 

1.5 

0.271*02 

0.271*53 

0.271*11 

0.27365 

0.25879 

0.25879 

2.0 

0.2671b 

0.26736 

0.26673 

0.26591 

0.25332 

0.25382 

2.5 

0.26009 

0.26010 

0.25958 

0.25863 

0.21*905 

0.21*905 

Odd  Mode 

\ e 

r\ 

0.0 

0.2 

Oj* 

0.6 

0.8 

1.0 

1.2 

0.0 

o.2d*56 

0.203011 

0.208682 

0.218826 

0.231*021 

0.253380 

0.281*1*62 

0.5 

0.23529 

0.  231*981 

0.239178 

0.21*9231 

0.266107 

0.286589 

0.317875 

1.0 

0.26968 

0.269917 

C. 271*278 

0.281*091 

0.29833 9 

0.316165 

0.31*2376 

1.5 

0.27U02 

0.271*031 

0.27891*0 

0.286029 

0.296391* 

0.308909 

0.327308 

2.0 

0.267m 

0.267602 

0. 271507 

0.276831 

0.283785 

0.292366 

0.301*981* 

2.5 

0.26009 

0.261693 

0.26JTO 

0.26781*1* 

0.27271*5 

0.27891*0 

0.28791*9 

CHAPTER  VH 


UfTEUPRETATIOK  OP  THE  CURVES  FOR  THE  QUALITY  FACTCF. 

(1)  Analysis  of  the  error  in  the  calculations. 

The  curve 3 far  tha  quality  factor  are  plotted  in  Figures  ? and 
6.  The  extremely  small  changes  in  the  quality  factor,  Q,  at  saall 
values  of  eccentricity  make  a careful  consideration  of  the  accuracy  range 
neoeasary.  The  limits  of  accuracy  in  the  determination  of  the  roots 
and  evaluation  of  the  integrals  was  noted  when  these  evaluations  were 
discussed  in  Chapter  VX.  It  is  now  desirable  to  see  hew  these  cocponent 
errors  are  reflected  in  the  final  result* 

The  error  calculations  are  simple  and  intermediate  steps 

are  omitted*  All  quotients  are  changed  to  product  form  by  expressing 

the  denominator  with  a negative  exponent.  AH  quantities  to 

m/n 

fractional  exponents  are  put  into  the  form  (a  + e;  and 

than  expanded  by  a binomial  expansion  to  two  terms.  All  numbers  are 

rounded  to  two  significant  figures. 

Tha  error  in  the  determination  of  the  derivative  root 
value  was  assumed  to  be  0.0000?.  This  may  also  be  considered  as  the 
error  appearing  in  ^^s•  It  develops  that  this  error  is  the  dominating 
one*  It  is  therefore  designated  by  the  symbol  e*  while  all  the  other 
errors  are  denoted  by  multiples  of  e where  e also  signifies  an  error 
of  0.0000?. 


) 


To  determine  the  error  in  kr  i 


iP^Jl  + ( £/s  4 e»)2  P2*2] 
4 1.3e«)  67 


The  error  in  evaluation  of  the  Integrals  I and  II  was  assumed 
to  be  0.00005,  and  the  error  in  the  evaluation  of  the  integrals  I* 
and  H*  was  considered  as  0*0005*  The  first  of  these  errors  is 
designated  by  e and  the  second  by  10  e.  The  error  in  the  evaluation 
of  V can  be  written  ast 


v«in>  4-  rn 

* (I*  e)(ll4  lOe)  4 (I*  ^-10e)(II  4 a) 

^ I II*  f I'  n 4 66» 

The  error  in  the  calculation  of  the  integral  HI  was  assumed 
to  be  0.00005  so  thati 

v/m  ^ ~~  •66?. 

3 4 e 

- 2/3  4 28e 

The  error  in  the  coefficient  of  Q,  which  is  independent  of  R,  can 
then  be  written* 


ooef  Q s 


3 r HI  J(r)  z 
_ 2/3  ♦ 28e 

3 (2  ♦ #')(•£  4-s) 

-2/9  4*  .le*  4 9 *k  e 


The  equations  for  estimating  error  arc  not  a developmental  step  and 
for  this  reason  they  are  not  umbered 


_ .uuyjcr-E  -wafrr  *»y-*~ 
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Similarly  it  was  found  that  far  the  P.^  coefficient: 


coef  R2*  ? + ^e'  -f  3e 


and  for  the  Ir  coefficient: 


coef  R" 


•5  4 17«  + 13e< 


For  R s 1 


Q */s  ~ .3  * l5.1e'  4 19e 
- .3  - .0017 


(for  both  e’  and.  e taken  as 
0.0000^) 


For  R s 2.5 


Q 


•3  + 6e»  + 6.3  e 

.3  4 .0007 


(for  both  ef  and  e taken  aa 

0.00005) 


It  is  noted  that  the  error  in  the  determination  of  the  root 

introduces  about  the  same  error  in  the  final  result  as  the  sum  of 

all  the  other  errors  in  the  determinations  of  the  integrals.  If  a 

more  conservative  estimate  of  the  error  in  the  determination  of 
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the  root  is  made  of,  say,  0.0002,  the  estimated  error  in  the 
determination  of  Q 5 /a  becomes  O.OOii. 

X2)  SygggEC  oWMOacfc  curves 

Ihen  a problem  is  arranged  for  maximum  facility  in  numerical 
oonputatic!?*  the  significant  components  may  be  mixed  so  that  it 


£2 


oeccnes  difficult  to  see  the  relations  of  tho  Individual  changes 
that  cochins  to  cavje  the  overall  result.  Inspection  of  the 
curves  far  Q in  figures  £ and  6 shews  that  the  quality  factor 
varies  very  slightly  for  values  of  eccentricity  less  than  0*U,  hut 
a cursory  examination  of  ‘the  data  in  Tables  I,  II#  Hi  shows  that 
these  slight  changes  have  not  been  the  result  of  snail  changes  in  the 
o opponents,  but#  rather#  of  ocopensating  changes  of  ouch  larger  magnitudes 
than  the  resulting  change  in  the  quality  factor#  Q* 

To  examine  then  qualitatively#  the  Individual  variations  in 
the  stared  energy#  the  power  loos  in  the  end  walls#  and  the  power  loss 
in  tho  sidewall  are  plotted  in  Figure  7 to  slide  rule  accuracy*  The 
characteristics  depend  on  the  ratio  of  the  average  diameter  to  the 
length  so  that  it  is  necessary  to  choose  a fixed  value  of  Rj  the 
convenient  value  of  R ■ 1 is  chosen*  The  curves  in  Figure  7 are 
based  on  unit  magnitude  at  zero  eccentricity,  and  they  provide  no 
information  concerning  the  relative  magnitude  of  the  changes  which 
occur* 


The  volume  of  the  cavity  is  directly  proportional  to  the 
product  of  the  major  and  minor  axes.  For  a given  perimeter,  the 
cross  section  area  is  a maxima  at  *ero  eccentricity,  arr  It  may  seem 
odd  that  the  stored  energy  increases  when  the  volume  is  decreased* 

At  values  of  eccentricity#  the  volume  decreases  very  little  as 

the  circle  is  deformed;  'a  plot  of  the  volume  on  a per  unit  basis  is 


- t fir 
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included  in  Figure  7 for  comparison  with  the  cnange  in  the  stored 
energy.  The  increase  in  stored  energy  can  be  attributed  to  the 
change  in  resonant  wavelength,  and,  perhaps,  a slightly  nore  efficient 
distribution  of  the  fields  in  the  cavity  with  the  change  in  eccentricity 
Although  this  change  in  distribution  would  be  difficult  to  demonstrate 
mathematically,  it  seems  to  have  some  intuitive  Justification.  The 
liathieu  function  at  small  eccentricities  approximates  a slightly 
distorted  Bessel  function.  The  lack  of  symsietry  of  the  dominating 
first  order  Bessel  function  in  the  range  from  ssero  to  the  first 
derivative  root  makes  it  seem  possible  that  a slightly  non-symmetric 
cylinder  might  use  such  a function  more  efficiently  than  the  perfect 
cylinder  could. 

At  an  eccentricity  of  0.5  in  the  odd  mode,  it  in  noted  that 
the  stored  energy  is  still  increasing.  Certainly  this  increase  could 
not  continue  indefinitely  and  must  change  to  a decrease  when  the 
volume  begins  to  decrease  rapidly# 

(3)  TKwm^nn  of  Results  and  Conclusions 

The  numerical  values  of  Q for  aero  eccentricity  agree 

numerically  with  those  plotted  by  -u^lgomery  for  the  circular 
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cylindrical  cavity. 

The  plot  of  the  quality  factor  in  Figure  5 for  the  odd  mode  shows 
that  the  quality  factor  may  decrease  in  value  slightly  for  very  small 


aiwamts  of  eccentricity.  Tills  decrease  is  of  the  same  order  as 
the  possible  error  in  calculations  and  should  not  be  given  undue 
consideration*  It  seems  reasonable  to  conclude  that  the  quality 
factor  remains  constant  for  values  of  eccentricity  less  than  0,25 
for  excitation,  in  either  the  erven  or  odd  node* 

Ordinarily,  a deformed  circular  cylinder  will  be  excited 
at  the  same  time  in  both  the  odd  and  even  modes.  If  the  excitation 
orientation  can  be  controlled  relative  to  the  deformation,  it 
•will,  be  preferable  to  excite  the  even  mode,  since  both  the  wavelength 
and  quality  factor  change  less  in  that  mode. 

It  is  unfortunate  that  other  modes  could  not  be  evaluated  by 
an  approximate  means j most  of  the  individual  quantities  involved 

7 **v 

in  the  expression  for  the  quality  factor  were  calculated  by  Tang 
but  calculations  made  using  his  values  have  a possible  error  range  that 
in  larger  than  the  magnitude  of  the  change  in  Q.  Calculations  would 
be  sinplified  a little  for  modes  for  which  the  newer  Tables  Relating  to 
Mathieu  Functions  could  be  used  since  some  of  the  integrals  are 
evaluated  directly  in  those  tables. 
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